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We have studied the field-orientational dependence of zero-energy density of states (FODOS) for 
a series of systems with different Fermi surface and superconducting gap structures. Instead of 
phenomenological Doppler-shift method, we use an approximate analytical solution of Eilenberger 
equation together with self-consistent determination of order parameter and a variational treatment 
of vortex lattice. First, we compare zero-energy density of states (ZEDOS) when a magnetic field is 
applied in the nodal direction (u no d e (0)) and in the antinodal direction {v a nti{0)), by taking account 
of the field-angle dependence of order parameter. As a result, we found that there exists a crossover 
magnetic field H* so that v a nu(Q) > &W>de(0) for H < H* , while v no de(0) > v a nu(Q) for H > H* , 
consistent with our previous analyses. Next, we showed that H* and the shape of FODOS are 
determined by contribution from the small part of Fermi surface where Fermi velocity is parallel 
to field-rotational plane. In particular, we found that H* is lowered and FODOS has broader 
minima, when a superconducting gap has point nodes, in contrast to the result of the Doppler-shift 
method. We also studied the effects of in-plane anisotropy of Fermi surface. We found that in-plane 
anisotropy of quasi-two dimensional Fermi surface sometimes becomes larger than the effects of 
Doppler-shift and can destroy the Doppler-shift predominant region. In particular, this tendency is 
strong in a multi-band system where superconducting coherence lengths are isotropic. Finally, we 
addressed the problem of cusp-like singularity in YNi2B2C and present a possible explanation of 
this phenomenon. 

PACS numbers: 



I. INTRODUCTION 

While it is accepted that electron correlation drives 
non-s-wave-pairing superconductivity, determination of 
detailed superconducting gap structure still remains a 
challenge for both theory and experiment. On a the- 
oretical side, analyses starting at microscopic Hamilto- 
nians have succeeded in identifying gap structures of 
several materials. On an experimental side, it is get- 
ting recognized that field-orientational dependence of 
specific heat or thermal conductivity carries a consid- 
erable information on gap structure. Since the low- 
energy density of states changes according to the an- 
gle between gap nodes and magnetic field, one can de- 
rive gap structure by tracing a change of thermodynamic 
quantities as sweeping magnetic field in several direc- 
tions. Actually, this method has been applied to a series 
of superconductors including S^RuO^i^Siii^, CeCoIng££, 
K-(ET) 2 Cu(NCS)2i, YNi 2 B2G&2, PrOs 4 Sbiz±a and 
UT^Abjii. However, their results left serious problems 
yet to be overcome. For example, in the case of Sr 2 Ru04, 
thermal conductivity (k) does not show any change with 
rotating magnetic field within the ab-plane, while specific 
heat (C) shows clear four-fold oscillation. As to CeCoIns, 
k and C show incompatible behavior. In YM2B2C, both 
k and C show a cusp-like singularity as sweeping mag- 
netic field within the superconducting plane. 

In the previous paper, we addressed the problem of 
Si'2Ru04 in terms of the fine-structure of density of 
states 12 . Whereas, the latter two problems arc left un- 
solved. While it is proposed that the cusp-like singularity 
in YN12B2C is caused by point node in gap structure^, 



their discussion is phenomenological in nature and unre- 
liable, as we will show. 

The confusion as to the experiment is probably at- 
tributed to a lack of firm theoretical basis. So far, the 
experimental data has been usually interpreted in terms 
of the "Doppler-shift criterion" . This criterion says larger 
density of states is obtained when a magnetic field is ap- 
plied in antinodal rather than nodal directions. However, 
this criterion is obtained on the basis of the phenomeno- 
logical Doppler-shift method that can be justified only 
at a limit of low-magnetic field. Hence, it is not clear 
whether this criterion is still valid in the experimentally 
accessible higher magnetic fields. 

In fact, in the previous paper, we showed that there oc- 
curs a crossover from Doppler-shift predominant region 
at low magnetic fields to core states predominant region 
at high magnetic fieldsi^. In the core state predominant 
region, larger density of states is obtained for field in 
nodal direction, in contrast to the "Doppler-shift crite- 
rion" . Hence, in order to derive a correct nodal structure 
from the field- rotational experiments, it is highly desir- 
able to establish a theoretical analysis applicable more 
widely. 

To this end, we present a more quantitative analyses 
of the crossover behavior by taking account of the field- 
angle dependence of order parameter. In particular, we 
are interested in how the Fermi surface or superconduct- 
ing gap structures affect the crossover behavior. 

Furthermore, we would like to investigate field- 
orientational dependence of zero-energy density of states 
(FODOS) in detail at the Doppler-shift predominant re- 
gion. It is interesting to clarify how information on nodal 
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structure is obtained from the shape of FODOS. 

The structure of this paper is as follows. In the next 
section, we summarize our method. In section IIIII we 
present a quantitative study of the crossover behavior 
for several models with different Fermi surface and nodal 
structures. In section llVl we examine the effects of in- 
plane anisotropy of Fermi surface to the crossover be- 
havior. In section we address the problem of cusp- 
singularity in YNI2B2C. And section IVII is devoted to 
conclusions. 



II. METHOD 

In this section, we introduce an approximate analyt- 
ical solution of the quasiclassical Eilenberger equation, 
which was originally proposed by Peschi^, and recently 
extended to unconventional superconductors by Dahm 
et al^. With this method, one can obtain exact results 
at H = H C 2, and quantitatively reliable results to much 
lower fields. All our analyses in this paper are based on 
this approximate analytical approach. 

We start with the quasiclassical Eilenberger 
equationiSiSS. This equation provides a convenient 
method to analyze inhomogeneous state of supercon- 
ductors by decoupling slowly varying order parameter 
from fast oscillating Fermi-particle degree of freedom. In 
this paper, we limit ourselves to the case of spin-singlet 
superconductors in clean limit. Under these conditions, 
the quasiclassical Eilenberger equations read 



v F ■ (V - i— A(r)) / + 2uj n f ~ 2iAg = 0, (1) 
c 



v F • (V + i— A(r))/t - 2w n /t + 2iA*g = 0, (2) 



with the normalization condition 



g 2 - fP = 1 



(3) 



Here, g = g(r,k F ,iu> n ), f = /(r,k F ,KJ„) and /+ = 
/^(r, k F , ioj n ) are normal and anomalous components of 
quasiclassical Green function, lj tl — (2n + 1)ttT is Mat- 
subara frequency, and A = A(r, k F ) denotes order pa- 
rameter. The order parameter A is determined in a self- 
consistent way, by solving the following gap equation. 

A(r,k F ) = -27r£7V n (0) ^ <V(k F ,k' F ) 

0<LU n <UJ c 

xf(r,k' F ,iu n ) >FS(k' F ) • (4) 



Free energy Q s can also be expressed with the quasi- 
classical Green function in a closed formi2i2i. 



fB — H) 2 
dr± — '- + mTu n (Q) 



8tt 



0<w„<oj c 



£ <5±l (A */ + A/t)> F5 , (5) 



with free energy in the normal state, H, external 
magnetic field, and B, total magnetic field. Equations 
form a basis to study inhomogeneous supercon- 
ductors. 

Next, we introduce an approximate analytical solution 
of Eilenberger equation. Equations Q-Q can be solved 
analytically only at H = H C 2- However, by extending the 
analytical solution to lower fields, one can obtain quan- 
titatively reliable results for density of states etc, well 
below H C 2- To obtain approximate analytical solution 
at arbitrary fields, we have to make three assumptions 
which are exactly justified at H = H C 2- 

First, we decouple the order parameter into the spatial 
and momentum parts 



A(r,k F ) = A *(r)$(k F ), 
and normalize ^(r) and <i>(k F ) as 
< |*(r)| 2 > r =l, 

< |$(k F )| 2 > FS = 1. 



(G) 



(7) 



(8) 



Here, < • ■ ■ > r means averaging in the real space. If 
\l/(r) is periodic, this is equivalent to taking average over 
a unit cell. Then, we assume the spatial variation ^(r) 
is described by Abrikosov vortex function, namely, 



*(r) = <Mr) = l J2 



N 



exp 



2nin 



y-^yi 



n— — 00 



17 5 



yo 



{x-nx ) 2 . (9) 



^a(^) can be obtained by solving linearized Ginzburg- 
Landau equation with different coherence lengths in x 
(£ x ) and in y directions. Here, we set x and y axes 
perpendicular to the magnetic field. N — (jf^)^ is a 
normalization constant to make < ^(r)! 2 > r = 1. Area 
of unit cell, A, can be related to the total magnetic field 
5 by 



7TC 

A = x y = —. 

eB 



(10) 



Here, u n (Q) = J dQ± is Zero-energy density of states 6 denotes the anisotropy of superconducting coherence 



FS 



(ZEDOS) in the normal state, w c , a cut-off energy, 

V(k F ,k' F ), a paring potential, and < ••■ > F s means 

taking average on the Fermi surface, namely, < • • • > F s= 
r dfl 1 

J !/„(0) Vp ' ' ' ' 



lengths, i.e., S — j L . One can estimate S to be S 



for single-band superconductors. 



FS 



<»f,> 

Here we note that \I>(r) denotes an Abrikosov vortex 
lattice which is spanned by two lattice vectors, (0, yo) and 
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(xo,yi). When yi — 0, ^(r) describes a square (rect- 
angular) lattice, while when j/i = \&(r) describes a 
triangular lattice. However, the shape of the unit cell is 
unimportant in the subsequent analyses. 

As a second assumption, we set Ginzburg-Landau pa- 
rameter k to be much larger than 1, hence B = H. This 
assumption is justified for a lot of unconventional super- 
conductors including YNi 2 B 2 C. 

The third assumption is that we ignore spatial fluctu- 
ation of g(r, kp, iu> n ), namely, < g 2 > r —< g > 2 . Since 
at H — H C 2, the order parameter vanishes and trans- 
lational symmetry is recovered, this assumption is also 
justified. The validity of this assumption has been justi- 
fied for lower fields^. 

With these three assumptions, one can obtain an ap- 
proximate analytical solution of eqs. Q~© a ^ arbitrary 
magnetic fields by utilizing the operator techniques intro- 
duced in ref. [z^. We summarize the basic results below. 
For details, see, for example, refs. ^| and fill 



[1 ] spatial average of g(r, k F ,iu/ n ) 
< g(r,k F ,iu} n ) > r = -Re 



1 



s/l + P(v F ,iu n ) 



(11) 



where, 

P(v^„) = ^|A |>(ffiM)* (1 



T( |V FJ J 

F(A,v F± ,u n )), (12) 



OO 

F(A,v F ±,cj n ) = J duexp(- 



Tr\VFA_\\ 2 



u 2 ). (13) 



Here, effective Fermi velocity vp± is expressed as 



= \f8v Fx e x H — 



(14) 



with the projection of Fermi velocity on x (v Fx ) 
and y (v Fy ) axes. 

[2 ] DOS v{e) 
Ke) 



,(0) 



= - « g(r,k F) iuj n ->• e + iO) > r > F s ■ (15) 



In particular, the formula of ZEDOS becomes quite 
simple. 



1/(0) 



MO) '^ 1 + M|Aop(J^fl) a 



>FS ■ (16) 



[3 ] gap equation 

In order to obtain gap equation, we decouple 
V(k F ,k' F ) as V(k F ,k' F ) = V<S>*{k F )$(k' F ), then 



we obtain 

F{A,v F± ,w n ) 
VI + P(-v F ,iuj n ) 

where 7 = 0.5772 • • • is Euler's constant. The for- 
mula to determine H C 2 can be obtained by setting 
A = in eq. l(T7|) . 



, 2eW = r 
irT r / sinh s 



l ~^^\ <m k F )\ 2 



[4 ] free energy 



! V F-L I n 2 , x 

Xe S7rT 2 A > FS . (18) 



2nT c 



n fl -n„ = -//„(o)|A | 2 V f^<|$(k F )| 2 



0<LO n <UJ c 

P(v F , iw„)F(A, v f _l, u;„) 



V / l + J P(v F ,iWn)(l + v / l + -P(v J F,^„)) 2 



>FS ■ (19) 



We calculate density of states by the following procedure. 
First, Fermi velocity v F and nodal structure ^ki?) are 
determined by assuming a phenomenological model or 
starting at a microscopic Hamiltonian. Then, at a given 
temperature T and a magnetic field H, one can calculate 
order parameter amplitude Ao by solving eq. (|17fl . By 
substituting A into eq. ijTlfl) . one can obtain DOS v{e). 
When magnetic field is applied off a symmetry axis, the 
ratio of coherence lengths 8 deviates from 1. In this case, 
we determine 8 variationally by minimizing free energy 
eq. i|19|) with respect to 8. This prescription is equiva- 
lent to the variational ansatz of Abrikosov vortex lattice 
introduced in refs. |16| and I23L 

One can incorporate multi-band properties into this 
approximate analytical approach with a slight modifica- 
tion of [l]-[4}iSiii. In multi-band systems, k F , \r F , /, f\ 
g, A, v n (0) become band-dependent. To specify band 
index, we write these quantities as k F , v F , f a , f a \ g a , 
A a , ^"(0). Next, we decouple order parameter as 



A Q (r,k F ) = A£* Q (r)$ Q (k F ). 



(20) 



Here, we assume ^"(r) to be described by Abrikosov lat- 
tice, ^aM- Furthermore, we assume that the anisotropy 
of coherence lengths S a are the same for all bands (8 a — 
8), and determine 8 by minimizing free energy. Due to 
inter-band mixing effect by a pair-hopping term, the spa- 
tial variation of A of each bands are not independent any 
longer. Then, one can take account of the inter-band 
mixing effect through minimization of 8. 

As we will discuss in section. we are interested in 
a system where quasi-two-dimensional Fermi surface and 
three-dimensional Fermi surface co-exist. In such a sys- 
tem, superconducting coherence length may be isotropic 
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due to the three-dimensional Fermi surface, and lead to 
isotropic vortex structure. In this case, one may ex- 
pect unusual phenomena characteristic of multi-band sys- 
tem, in the sense that quasi-two-dimensional Fermi sur- 
face and isotropic vortex structure are incompatible in a 
single-band superconductor. 

Here, we note that some complication may arise in 
vortex structure due to multi-band effects. For example, 
it is proposed that order parameter of different bands 
may have zero-points at different places. However, we 
believe our results will not be influenced by the detail of 
vortex structure and ignore them. 

Here, let us check how [l]-[4] arc modified in a multi- 
band case, [ljspatial average of g(r, kp , iuj n ) also holds 
for multi-band systems, by making each quantity in eqs. 
(HTJ-CU band-dependent. [2]DOS v(e) and [4]free en- 
ergy are also correct by making the quantities band- 
dependent, and summing over band indices in the right 
sides of eqs. (TToT) . and lfT§)l. Only gap equa- 

tion in [3] needs substantial modification. In a multi- 
band system, we decouple pair potential V a/3 (kp, k'p) as 
V{k F ,k^) = V al3 ^*(k F )^(kf), then gap equation eq. 
(If 7|l becomes 

A «log(^)=£^(0)V^ £ ^ 

0<UJ„<UJ c 

x < |**(kg)| a f ^ 



*(k F ) = y— (1-sin 2 6» cos Acfi) . (27) 

These models are often used to approximate more com- 
plicated realistic materials. Hence, it is important to 
study the basic properties of these models in detail. 

Throughout our paper, we specify the direction of mag- 
netic field with two parameters 9q and (fio as 

H = H(s'm9o cos(/>oe Q + sin^o sin^o e h + cos 8oe c ). (28) 

In principle, we fix 8q = 90° and study the change of 
ZEDOS with varying (fio in this paper. In other words, 
we consider the ab-plane as a field-rotational plane. As to 
the three models above, (fio = 0° and (fio = 45° correspond 
to the cases where a magnetic field is applied in nodal 
direction and antinodal direction, respectively. 

Here, we focus on ZEDOS at low temperatures, be- 
cause universal behaviors are expected in this region. We 
note that some cares are necessary to discuss thermody- 
namic quantities at higher temperature, because DOS at 
high energies depend on the magnetic field direction quite 
differently from ZEDOS±2. Thus we fix T = 0.1T C and 
lu c = 20irT c throughout this section. At this tempera- 
ture, ZEDOS well represents the thermodynamic prop- 
erties of a system. 



A. Crossover behavior of ZEDOS and order 
parameter 



III. QUANTITATIVE STUDY OF THE 
CROSSOVER BEHAVIOR 

In this section, we present a quantitative study of the 
crossover behavior for the following three models. Below, 
we take e a , and e c as a right-handed orthonormal 
basis. 

(a) line nodes on spherical Fermi surface 
Vf = vp (sin 9 cos <fie a + sin 9 sin (fief, + cos9e c ). (22) 



$(k F ) = y -j sin 2 (9sin20. (23) 

(b) line nodes on cylindrical Fermi surface 

Vp = vf(cos <fie a + sin^e?, + esirik z e c ). (24) 

$(k F ) = V2sin20. (25) 

Here, e gives a small c-axis dispersion. We set e = 
0.05. 

(c) point nodes on spherical Fermi surface 

•vp = vf (sin 9 cos (fie a + sin 9 sin (fie^ + cos9e c ). (26) 



First, we consider ZEDOS when magnetic field is ap- 
plied in the nodal direction {v no de{fy) and in the antin- 
odal direction (v ant i(Q)) . In the previous paper—, we 
showed that the Doppler-shift criterion is correct at low 
fields, while it is broken at high fields due to the quasi- 
particles propagating parallel to magnetic field. In the 
previous paper, we used an assumption that Ao does 
not depend on field direction. In the present paper, we 
calculate ZEDOS with determining Ao self-consistcntly. 
In general, Ao changes as field direction varies. In par- 
ticular, near H C 2, in-plane anisotropy of H C 2 influences 
FODOS through the field-angle dependence of A . Per- 
haps, the anisotropy of H C 2 continues to affect FODOS 
to lower fields, and suppress the Doppler-shift predom- 
inant region. Hence, it is of interest to estimate how 
the field-angle dependence of Ao modifies the crossover 
behavior. 

In fig. ^ we show the FODOS for the three models 
introduced above. These show a similar crossover be- 
havior with different crossover fields H* . For H < H* , 
v a nti(0) > ^node(O), consistent with the "Doppler-shift 
criterion", while for H > H* , v no de(0) > v an u(0)- These 
results are consistent with our previous analyses^. Since 
H C 2 is smaller when field is in the nodal direction (H || n) 
than in the antinodal direction (H \\ a) in all the three 
cases, the anisotropy of H C 2 generally supports the rela- 
tion u node (0) > VantiiP) at high fields. 
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FIG. 1: Field dependence of ZEDOS under magnetic field in 
the nodal and antinodal directions for (a)line nodes on spher- 
ical Fermi surface, (b)line nodes on cylindrical Fermi surface, 
and (c)point nodes on spherical Fermi surface. (Inset) Field 
dependence of Ao. Ao is normalized with T c , and magnetic 
field is normalized with H C 2 when a field is in the antinodal 
direction. The agreement between (a) and the result of Eilen- 
berger approach (fig. 4 of ref. 24) is quite good. 
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FIG. 2: Field dependence of <5Ao and v a nti(Q) — v n ode(0) 
for (a)line nodes on spherical Fermi surface, (b)line nodes 
on cylindrical Fermi surface, and (c)point nodes on spheri- 
cal Fermi surface. <5Ao is normalized with T c , and magnetic 
field is normalized with H C 2 when a field is in the antinodal 
direction. 
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In order to examine how the anisotropy of Ao affects 
FODOS at low fields, we plot SA = A (H \\ a)-A (H \\ 
n) in fig. together with the anisotropy of ZEDOS, 
i / onti(0)— i / node(0). One can find that 8A$ changes its sign 
near H* and becomes negative for H < H* . This means 
that anisotropy of Ao has the same sign with the relation 
fanti(O) > Vnode(O) in the Doppler-shift predominant re- 
gion. Combining this with the effects of anisotropic H C 2 
at high magnetic fields, it is shown that our previous 
result on the crossover behavior is not altered by tak- 
ing account of the field- angle dependence of A . 

Next, let us study individual cases. First, let us com- 
pare the cases of line nodes on spherical Fermi surface(fig. 
EJa)) and on cylindrical Fermi surface(fig. HJb)). In spite 
of the difference in dimensionality of Fermi surfaces, both 
cases show similar behavior of v no de{0) and f an ti(0) with 
very close H* . In fact, the small c-axis dispersion of the 
cylindrical Fermi surface plays a key role in this simi- 
larity. In eq. Q16|). information of Fermi surface is in- 
cluded only in |vfj_|. In the case of cylindrical Fermi 
surface, |vfj_| 2 oc 5cos 2 (0 — (f>o) + |e 2 sin 2 k z . While 
e is much smaller than 1, one cannot neglect the term 
^e 2 sin 2 k z . In the case of anisotropic Fermi surface, one 



has S 



<vj> 



e in the whole range of mag- 



£« V <<*> 

netic field, as a result of minimizing the free energy eq. 

Hence, the in-plane and the c-axis components of 
effective Fermi velocity, <5cos 2 (</>— </>o) and |e 2 sin 2 k z are 
scaled to be equivalent in magnitude. 

Next, we study the case of point nodes on spherical 
Fermi surface (fig. Efc)). In this case, one can find that 
the crossover field H* is smaller (H* ~ Q.2H C 2) than 
the case of line nodes {H* ~ OAH^). This difference 
in H* can be understood in terms of the shape of gap 
minima. To see this, we note that the reversal of the 
Doppler-shift criterion at high magnetic fields is caused 
by the quasiparticles propagating parallel to magnetic 
fielcU^. At high magnetic fields, contribution from such 
quasiparticles to ZEDOS, vqp\\h, is smaller for H \\ a 
than H || n, resulting in the relation v n ode(0) > v a nti{fy- 
Here, for the point node case, |$(kp)| is quadratic near 
gap minima (for example, at 9 = 90°, <j> = 0°), while 
it is linear in the line node cases. Therefore, the order 
parameter is smaller near gap minima in the point node 
case, leading to larger vqp\\h f° r H || n. On the other 
hand, there is not much difference in vq P \\ h when H \\ a. 
As a result, v n ode{fy — ^anti(O) is larger in the point node 
case and continues to be positive to lower field. 

One may suppose that the smaller H* can be at- 
tributed to the gap structure in the "polar part" (9 ~ 0°). 
Since in the point node case, |$(kjr)| ^ except for 
the "equatorial part" (9 ~ 90°) of the Fermi surface, 
the effect of Doppler-shift is weaker in the "polar part" , 
and may lead to the smaller H*. However, this does 
not hold. In fact, if we replace the gap function eq. 
(|T7Jl with <f>(k F ) = yl — sin 2 9 cos 40, i.e., a point node 
model which varies linearly near gap minima, we find 
that H* ~ 0.4iJ C 2(fig. EJ, comparable to the line node 
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FIG. 3: Field dependence of SAo and v ant i(0) — ^node(O) for 
spherical Fermi surface with point nodes. Gap function is 
given by $(k_F) = yl - sin 2 8 cos 4<f>. 



cases, although |$(kp)| is larger than the quadratic point 
node model (eq. everywhere on the Fermi surface. 



B. Field-orientational dependence of ZEDOS 

In this subsection, we study not only the amplitude 
but also the whole field-orientational dependence of ZE- 
DOS. In particular, the shape of minima is of our in- 
terest, since we will discuss the origin of cusp-like min- 
ima which has been attributed to a point node in the 
Doppler-shift methodic. In fig. 0] we show our results 
for the three models. These calculations are done at 
H = 0.lH c2 . One can find similar FODOS for the line 
node cases (fig. IHa),(b)) with a slight difference in am- 
plitudes. Whereas in the point node case (fig. EJc)), 
one obtains FODOS with broader minima than the line 
node cases. Furthermore, in the point node case, tiny 
maxima appear when magnetic field is applied right in 
the nodal direction (<fr = 0°). These features again re- 
flect the quadratic variation of order parameter near the 
gap minima. For comparison, we show FODOS for the 
gap function ^(ki?) = yT — sin 2 9 cos 40 in fig. where 
one can find narrower minima similar to the case of line 
nodes. In both cases, our calculation does not support 
the results of Doppler-shift method. 

For the better understanding of the features in FO- 
DOS, we plot contribution to FODOS from several parts 
of the Fermi surface in fig. From these figures, one can 
find that FODOS is determined by the contribution from 
8 r~j 90° on the spherical Fermi surface or k z ~ 0° on the 
cylindrical Fermi surface. In other words, Most contribu- 
tion comes from the part where effective Fermi velocity 
v_F_i_(eq. (|f 4|l ) is parallel to the field-rotational plane. 
One can coherently understand the features of FODOS 
and H* obtained in this section from this viewpoint. 

First, the difference between point node and line node 



7 




FIG. 4: Field-orientational dependence of ZEDOS for (a)line 
nodes on spherical Fermi surface, (b)line nodes on cylindrical 
Fermi surface, and (c)point nodes on spherical Fermi surface. 
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is unimportant, because these gap structures behave sim- 
ilarly at 9 ~ 90°, where Vp± is parallel to the field- 
rotational plane. The dimensionality of Fermi surface is 
also unimportant. While Fermi velocity is almost parallel 
to the field-rotational plane on a cylindrical Fermi sur- 
face, the effective Fermi velocity is so only near k z ~ 0°, 
due to small <5(~ e). Whereas it is important whether 
order parameter varies linearly or quadratically near gap 
minima, because this difference lies in 8 ~ 90°. In this 
sense, H* and FODOS are determined by the local struc- 
ture of order parameter on the part where vp± is parallel 
to the field-rotational plane, rather than the global nodal 
structure of the entire Fermi surface. 



IV. EFFECTS OF IN-PLANE ANISOTROPY OF 
FERMI SURFACE 

In this section, we examine the effects of in-plane 
anisotropy of Fermi surface on FODOS. The anisotropy 
of Fermi surface affects ZEDOS through \vpj_\ in eq. 

If the distribution of |v_pj_| changes according to 
the magnetic field direction, ZEDOS also depends on the 
field direction, even without gap nodes. In this section, 
we focus on the shift of H* with increasing anisotropy of 
Fermi surface. 

In order to take account of the in-plane anisotropy of 
Fermi surface, we consider a simplified model shown in 
fig. [3 We make Fermi surface by connecting the sets of 
wave numbers, i.e. by adding nested parts to a cylindrical 
Fermi surface. 
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(k a ,h) = (±(u + v),k b ), (-v<k b <v), (29) 



(k a ,k b ) = (k a ,±(u + v)), {-v<k a <v), (30) 



(k a , kf,) — (±v ± ucos (j>, ±v ± usincj)), (0 < <j> < ^)(31) 

Here, £ = v/u is a measure of anisotropy. When £ = 0.5 
and 1.0, the linear parts occupy 24.1 and 38.9% of the 
Fermi surface, respectively. 

We assume the absolute value of in-plane components 
of Fermi velocity vf are the same all over the Fermi sur- 
face, and set a c-axis component vfc — 0-05vf sin k z . We 
assume a d-^^^-like pairing, i.e. 



QQlf) cxcos(2tan- 1 (^)), 



(32) 



i.e., gap nodes are on the cylindrical parts of the Fermi 
surface. 

When magnetic field is applied perpendicular to one of 
the nested parts, \vf±\ ~0 on the nested parts, leading 
to suppression of ZEDOS. Whereas, when magnetic field 
is applied in the nodal direction, the nested parts does 
not influence ZEDOS so much. Hence, it is expected that 
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FIG. 6: Contribution to FODOS from several parts of the 
Fermi surface for (a)line nodes on spherical Fermi surface, 
(b)line nodes on cylindrical Fermi surface, and(c)point nodes 
on spherical Fermi surface. In these three cases, total FODOS 
is given by averaging over 9 (0° < 6 < 90°) or k z (~n < k z < 

7T). 
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FIG. 8: Comparison of u no d e (Q) and v a nti(0) for £ = 
FIG. 7: A schematic picture of the Fermi surface with in- 0,0.5, 1.0. Magnetic field is normalized with H c2 when afield 
plane anisotropy seen from the c-axis. Position of gap node i s m t ne antinodal direction, 
is shown with filled circles 

the Doppler-shift predominant region is suppressed as £ 
is increased. In fig. [HJ we show v an ti(0) — v n ode(0) for 
several values of C- One can find that H* is lowered as £ 
is increased. At £ = 1.0, the Doppler-shift predominant 
region is suppressed to H ~ O.O5_ff C 2(0o = 45°), i.e., the 
almost entire field range at which experiments have been 
done. 

This result suggests that many cares should be taken 
to apply the "Doppler-shift criterion" when Fermi sur- 
face is highly anisotropic. At least one should check that 
maxima and minima of FODOS are exchanged at high 
fields and low holds, in order to confirm that the system 
enters the Doppler-shift predominant region at low mag- 
netic fields. It seems that recent experiments do not find 
any exchange of maxima and minima^. 



V. ORIGIN OF THE CUSP-LIKE 
SINGULARITY 



In this section, we study origin of the cusp-like sin- 
gularity in FODOS observed through specific heal£ and 
thermal conductivity in YM2B2C. In both cases, the 
cusp-like minima are observed when magnetic field is ap- 
plied parallel to a or b-axis. So far, the cusp-like min- 
ima are attributed to point nodes on the superconducting 
gap. In particular, s+g-wave gap structure has been pro- 
posed for YNi2B2C 13 . However, as we showed in section 
IIIII rather broader minima of FODOS are obtained from 
s+g-wave-like quadratic point nodes. 

Here, we would like to propose a new mechanism that 
leads to cusp-like minima in FODOS. In our theory, two 
conditions are required for the appearance of cusp-like 
minima. The first condition is that a system has to have 
a quasi-two-dimensional nested Fermi surface. And the 
second condition is that the superconducting coherence 
lengths are isotropic. Under these conditions, cusp-like 
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minima appear when a magnetic field is applied perpen- 
dicular to the nested part of the Fermi surface. 

Here, we briefly describe how cusp-like minima appear 
from the above conditions, in terms of eq. I|16|) . Provided 
that superconducting coherence lengths are isotropic, one 
obtains S = ~ 1. When the magnetic field is applied 
at an angle <po with the normal line of the nested part of 
the Fermi surface, we have 



I v f_l| 



e 2 sin 2 k z ). 



(33) 



on the nested part, with a c-axis component of Fermi 
velocity, e sin k z . If the c-axis dispersion is weak enough 
(e << 1), one can ignore the second term in the right- 
hand side, and the contribution to ZEDOS from the 
nested part of the Fermi surface v nest (0) is written as 



V nest (0) OC 



2A| Ao |2(M^ 

7T ' u 1 \Vp Sill 0o ' 



2A|A $(k F )|- 



(34) 



for small aV Hence, FODOS has cusp- like minima at 
(fro = 0, i.e., when a magnetic field is applied perpendic- 
ular to the nested parts. 

The | sin <fio | dependence in eq. (|34(l allows one a 
simple physical interpretation. Suppose a quasiclassi- 
cal trajectorji22iSiSi of nested Fermi surface, which is 
running in the direction at an angle 0o with magnetic 
field, and passing near vortex cores (fig. EJ- On this 
trajectory, phase of order parameter changes drastically 
near vortex cores, leading to formation of Andreev bound 
states. Since the vortex cores are separated by the dis- 
tance r~ 



VA 



I gin ^ | on the trajectory, the density of Andreev 

bound states is proportional to ' sl ^° , hence make a 
cusp-like singularity in FODOS. This brief discussion re- 
veals that the cusp-like singularity is caused by the core 
states, which are not properly taken into account in the 
Doppler-shift method. 

Here, we discuss the basic properties of YNi2B2C. 
YNi 2 B 2 C has three bands (17th,18th,19th) crossing 
Fermi energji2^2SiSi. The 17th band, which mainly con- 
sists of Ni(4d) orbits, forms an open Fermi surface in 
the c-axis direction, hence has a quasi-two-dimensional 
na t ure 25 4 26 i 27 jt, j s re p 0r ted that the corresponding 
Fermi surface of LuNi 2 B 2 C has a nested part, with nest- 
ing vectors q ~ (0.56 x 27r)e a ,(0.56 x 27r)ef, 2 I. Kohn 
anomaly in phonon spectrum22*22iSi provides another ev- 
idence of the nesting. On the other hand, the 18th and 
the 19th bands form closed Fermi surface in the c-axis 
direction, hence have three-dimensional tendency. As to 
the gap structure, existence of nodes have been confirmed 
by specific heali^i and ultra-sonic attenuation^ etc. 

In order to discuss the features of YNi 2 B 2 C, we study 
the following two-band model. Our model consists of 
(a) a partly nested quasi-two-dimensional Fermi sur- 
face as shown in fig. [7J and (/3) a three-dimensional 



00 




FIG. 9: (upper)A schematic picture of a trajectory running 
at an angle <f>o with magnetic field. Vortex cores are shown 
with dashed lines, (lower) Andreev bound states formed on 
the trajectory. 
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FIG. 10: (a)FODOS of the two-band model. Gap function 
(i) is adopted for the nodal structure of the spherical Fermi 
surface, (b) Contribution to FODOS from each Fermi surface. 



FIG. 11: (a)FODOS of the two-band model. Gap function 
(ii) is adopted for the nodal structure of the spherical Fermi 
surface, (b) Contribution to FODOS from each Fermi surface. 



spherical Fermi surface. We suppose the normal-state 
density of states of these two Fermi surfaces are the 
same. We set ( — v/u — 0.137 for the first Fermi sur- 
face. Then, about 4% of the entire Fermi surfaces are 
nested, in agreement with ref. |23, where the propor- 
tion of the nested part is reported to be 4.4 ± 0.5%. As 
to the gap structure, we assume isotropic gap for the 
quasi- two-dimensional Fermi surface, and the following 
two cases for the spherical Fermi surface, with differ- 
ent nodal position: (i)<I>(ki?) = ^/^sin 2 9sin2(f> and (ii) 

$(ki?) = sin 2 6 cos 2</>. As to the interaction poten- 
tial V af3 in eq. we assume V af3 /V aa = 0.5 and 
we assume spherical Fermi sur- 
face has the larger superconducting instability, in order 
to make superconducting coherence lengths of this sys- 
tem isotropic^. 

In fig. filial and fig. fTT^ a). we give FODOS of this 
model. Calculation is done at H = 0.lH C 2((j>() = 0°). 
One can find a clear cusp-like minima in both cases. Fur- 
thermore, it is worth noting that FODOS has minima for 
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FIG. 12: Contribution to FODOS from several k z of the 
quasi-two-dimensional Fermi surface. This figure should be 
compared with fig. H^b) 
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H || a(<fio = 0), irrespective of the nodal position. In fig. 
ITUTb) and fig. [TTTbh we show contribution to FODOS 
from each Fermi surface. One can find that the FODOS 
arising from the quasi-two-dimensional Fermi surface is 
larger and dominate the contribution from the spherical 
Fermi surface. It is surprising that only tiny anisotropy 
of Fermi surface destroys the Doppler-shift predominant 
region below H = O.Lff^^o = 0°). It makes a contrast 
with the single band case in section llVl where as much 
as 40% of Fermi surface has to be nested to suppress the 
Doppler-shift predominant region to H ~ 0.05H c2 . 

One can understand the origin of the enhancement of 
the contribution from the quasi-two-dimensional Fermi 
surface, on the basis of the scenario in section II I II In 
section Mil we noted that largest contribution to FO- 
DOS comes from the part of Fermi surface where effec- 
tive Fermi velocity vp± is nearly parallel to the field- 
rotational plane. In the case of cylindrical Fermi surface 
in section Mil such parts are limited to k z ~ due to 
anisotropic coherence lengths 6 ~ e. However, in the 
two-band model in this section, the superconducting co- 
herence lengths are isotropic (<5 ~ 1), hence vp± is al- 
most parallel to the field-rotational plane all over the 
quasi-two-dimensional Fermi surface. Therefore, contri- 
bution to FODOS come from all k z and becomes much 
larger than the contribution from the nodal structure on 
the spherical Fermi surface. We plot contribution to FO- 
DOS from several k z in fig. [21 In this figure, one can 
find that all k z contribute to FODOS equally. 

The cusp-like singularity shown here is basically ow- 
ing to the co-existence of two-dimensional and three- 
dimensional Fermi surfaces and to the nesting nature 
of two-dimensional Fermi surface. This example shows 
that many cares are needed for an analysis of experimen- 
tal data in the multi-band system. Then, anisotropy of 
quasi-two-dimensional Fermi surface is enhanced and in- 
validate the simple argument based on the Doppler-shift 
methods. 



VI. CONCLUSIONS 

We study the influence of superconducting gap and 
Fermi surface structures on the ficld-oricntational depen- 



dence of ZEDOS. First, we studied the crossover behavior 
of ZEDOS, by taking account of the four-fold oscillation 
of A , and found that there exists a crossover magnetic 
field H* so that v anti (ti) > v node (Q) for H < H* , while 
v n ode{fy > Vanti{fy for H > H * , consistent with our pre- 
vious analyses^. 

Next, we have investigated the roles of gap structure 
and nature of Fermi surface. The crossover behavior and 
FODOS at the Doppler-shift predominant region have 
been discussed. We found that there are no significant 
differences in neither the value of H * nor the shape of FO- 
DOS between the spherical Fermi surface and the cylin- 
drical Fermi surface, in spite of the difference in the di- 
mensionality of Fermi surface. We also found that when 
quadratic point nodes exist, H* is lowered than the line 
node cases and FODOS has broader minima. These re- 
sults can be coherently understood by recognizing that 
contribution to FODOS mostly comes from the part on 
the Fermi surface where effective Fermi velocity is nearly 
parallel to the field-rotational plane. 

Next, we examine the effects of in-plane anisotropy of 
Fermi surface. We showed that the Doppler-shift pre- 
dominant region may be suppressed to lower fields by 
the anisotropy of Fermi surface. This tendency becomes 
especially strong in a multi-band superconductor which 
has nearly isotropic coherence lengths. Therefore, some 
cares are necessary for analyzing the field-rotational ex- 
periments for YM2B2C and CeCoIns. On the other 
hand, the field-rotational experiment and its interpre- 
tation based on the Doppler-shift criterion are powerful 
tool for such material as Sr 2 Ru04, Nao.35CoO2l.3H2O 3 ' 
and organic superconductors, which are well described 
by a single Fermi surface or have plural Fermi surfaces 
with the same dimensionality. 

Finally, we propose a novel mechanism of the cusp-like 
minima found in YM2B2C, in terms of Andreev bound 
states. The cusp-like minima are attributed to nesting in 
quasi-two-dimensional Fermi surface and the isotropy of 
superconducting coherence lengths. 

We are grateful to Y. Matsuda, T. Sakakibara, Y. 
Kato, K. Izawa, H. Kusunose and T. Watanabe for fruit- 
ful discussions. We would also like to thank K. Kamata 
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periments of YM2B2C. 
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